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ABSTRACT 

We compare model results from our semi-analytic merger tree based framework for 
high-redshift (z ~ 5 — 20) galaxy formation against reionization indicators including 
the Planck electron scattering optical depth {ves) and the ionizing photon emissivity 
(Pion) to constrain the particle mass of Warm Dark Matter (WDM). Our framework 
traces the Dark Matter (DM) and baryonic assembly of galaxies in 4 DM cosmolo¬ 
gies: Cold Dark Matter (CDM) and WDM with a particle mass of rux = 2.25, 3 and 
5 keV. It includes all the key processes of star formation, supernova feedback, the 
merger/accretion/ejection driven evolution of gas and stellar mass, and the effect of 
the ultra-violet background (UVB) created during reionization in photo-evaporating 
the gas content of galaxies in halos with My fs, lO^M©. We show that current Planck 
values rule out ^ 2.5 keV WDM, even in the physically unlikely scenario that 
all ionizing photons produced by these galaxies escape and contribute to reionization 
(i.e. /esc = !)■ With the largest number of UVB-suppressed galaxies, CDM faces a 
“stalling” of the reionization process with this effect decreasing with the disappear¬ 
ance of small-scale structure with decreasing rux- Finally, we find the bulk of the 
reionization photons come from galaxies with a halo mass My iS IO^Mq, stellar mass 
M* ;$ IO^Mq and UV magnitude —18 ;S Mjjy ^ — 13 in CDM. The progressive sup¬ 
pression of low-mass halos with decreasing rrix leads to a shift in the “reionization” 
population to larger (halo and stellar) masses of My ^ 10®Mq and ^ IO^Mq for 
irix ^ 3 keV WDM, although the UV limits effectively remain unchanged. 

Key words: Cosmology:dark matter-reionization-cosmological parameters-cosmic 
microwave background; galaxies-intergalactic medium 


1 INTRODUCTION 


According to the standard Lambda Cold Dark Matter 
(ACDM) cosmological model, galaxy formation proceeds hi¬ 
erarchically through time, driven by a cold (dark) matter 
of unknown composition. In this era of precision cosmology, 
its energy density has been measured to have a value of 
= 0.133 with baryons comprising ~ 0.0196 of 
the total ( Planck Collaboration et al.|2014'b K CDM clusters 
on all scales and has been remarkably successful in predict¬ 
ing the large scale structure of the Universe, the tempera¬ 
ture anisotropies measured by the Cosmic Microwave Back¬ 


ground (CMB) and Lyman-a forest statistics (e.g. Peebles 
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linshaw et al.||2013 Planck 


recently reviewed by Weinberg et al. (20131, CDM exhibits 


a number of small scales problems: it produces halo pro¬ 
files that are cuspy as opposed to the observationally pre¬ 
ferred constant density cores (Navarro et al. 1997[ [Subra- 
manian et al.|20OT l, it over-predicts the number of satellite 
and field galaxies as compared to observations (the “missing 
satellite problem”; Klypin et al. 1999 Moore et al.] 19991, 


predicts massive (Large Magellanic Cloud mass), concen¬ 
trated Galactic subhalos inconsistent with observations (e.g. 
Boylan-Kolchin et al.|20l2 1 and faces difficulty in producing 
typical disks due to ongoing mergers down to z ~ 1 (Wyse 


20011. The limited success of baryonic feedback in solving 
these small scale problems (e.g. Boylan-Kolchin et al.|20T2 
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Teyssier et al.|[2013 1 has prompted questions regarding the 
validity of the CDM scenario. A popular solution to the 
small scale problems in CDM cosmology involves invoking 
(~ keV) Warm Dark Matter (WDM) particles that erase 
small-scale power (Blumenthal et al.||1984 Bode et al.|2001 


e.g)Q Beyond its Standard Model, particle physics provides 
compelling motivation for WDM candidates, which can be 
as light as ~ 0(1) KeV, such as sterile neutrinos. These ap¬ 
pear in an extension of the Standard Model with 3 sterile 
neutrinos, out of which 2 could be as heavy as 1 — 10 GeV, 
while the lightest one could be within 0(1) keV (for a re- 
Abazajian et al.|2012|. There are some tantalising 


observational hints such as the 3.5 keV monochromatic line 
observed by XMM-Newton coming from the Perseus galaxy 
cluster. This line might arise from a light sterile neutrino an¬ 
nihilating into photons (Bulbul et al.||2014 Boyarsky et al. 
20141). 


Astrophysical constraints on the WDM particle mass 
range from rux Ss 3.3 keV using the Lya forest power spec¬ 


trum (Viel et al. 20131, rrix > 1.6 — 1.8 keV using num- 


ber c ounts of high-z gamma ray bursts (de Souza et al. 
2013), rrix ~ 1.3 keV using abundance matching of theo¬ 


retical and observed high -2 galaxies (Schultz et al. 2014[), 
rUx ~ 1 keV using dwarf spheroidal galaxy observations | de 


Vega & Sanchez 2010), simultaneously reproducing stellar 


mass functions and the Tully-Fisher relation for 2 = 0 — 3.5 


galaxies (Kang et al. 20131 and comparing the observed 


number density of 2 ~ 10 galaxies to tha t expected from the 
halo mass function (Pacucci et al.|2013l, and rrix ~ 0.5 keV 


inferred using the presence of supermassive black holes at 
2 ~ 5.8 (Barkana et al.|2001|. 


A combination of ground and space-based observa¬ 
tions has allowed a statistically signihcant data available for 


2 ~ 6 — 10 Lyman Break galaxies (LBGs; Oesch et al.|2010 

Bouwens et al.||20101 20111 Castellano et al.j 

2010 

McLure 

et al. 

2010| [Bradley et al. [2012[ [Oesch et al.| 

2013 

McLure 

et al. 

2013[ [Bowler et al.|2014b [Bouwens et al.|201'! 

Bowler 

et al. 

2014a 1 enabling such constraints to be extended to 


the high -2 Universe. Indeed, in Dayal et al. (2014bI, we pre¬ 


sented a merger-tree based semi-analytic model that traces 
both the DM and baryonic assembly of high -2 (2 ~ 7 — 15) 
galaxies in 4 cosmologies: CDM and WDM with particle 
masses m^, = 1.5,3 and 5keV. While the halo mass func¬ 
tion is indistinguishable for CDM and ^ 3 keV WDM 
down to halo masses as low as Mh — lO^ ^M©, small mass 
(Mh — 10®'®M©) structures are suppressed more severely in 
the 1.5 keV scenario. We showed that the ultraviolet lumi¬ 
nosity functions (UV LFs), mass to light (M/L ratios) and 
stellar mass densities (SMD) are the same in all the 4 models 
for the massive, luminous (Muv ^ —18) galaxies that have 
been observed so far. However, the delay in structure for¬ 
mation in the 1.5 keV WDM scenario leads to a delayed and 
accelerated baryonic assembly resulting in a steeper SMD 


^ However, some works caution that the extremely low mass 
WDM particles required to make constant density cores prevent 
the very formation of dwarf galaxies ( [Maccio et al.|2012| |Schnei-| 
der et al.|2014 1 . 


evolution with redshift, that should be detectable by the 
James Webh Space Telescope (J WST) integrating down to 
Muv — —16.5 at 2 10 (see also Calura et al.||20l4 

ernato et al.||2014l. 


Gov- 


Since low-mass high -2 galaxies are now believed to be 
the main sources of reionization photons in the early Uni- 


verse (e.g. Barkana & Loeb 

2001 Ciardi & Ferrara 2005 

Choudhury & Ferrara 2007 

Choudhury et al. 2008), their 


delayed and accelerated assembly would naturally lead to a 
corresponding delay and acceleration in the reionization his¬ 
tory. Now that we have verified that our theoretical galaxy 
populations are in agreement with observations at 2 ~ 7—12 
for rrix = 1.5, 3 and 5keV, we build reionization histories for 
4 DM models with = 2.25, 3 and 5keV to see if we can 
obtain tighter constraints on by comparing to existing 
observations of the Planck CMB electron scattering optical 
depth and the ionizing photon emissivity. In addition to in¬ 


ternal feedback from supernovae (see Sec. 2.3 Dayal et al. 


2014b I our model self-consistently includes the “external” 


feedback effect from reionization in photo-evaporating gas 


from DM halos (Sec. 2.2 1 , and yields both the reionization 


history, and the ionizing photon contribution across differ¬ 
ent halo mass, stellar mass and magnitude bins (Sec. [^, as 
explained in what follows. 

The cosmological parameters used in this 
work correspond to (Dm, Da, Db, h, n-s, ug) = 
(0.2725,0.702,0.04,0.7,0.96,0.83), consistent with the 


latest results from the Planck collaboration ( Planck Col 
[laboration et aD] |2014b[ ) and we quote all quantities in 
comoving units unless stated otherwise. 


2 THEORETICAL MODEL 


We now briefly summarise the theoretical model and inter¬ 
ested readers are refereed to Dayal et al. ( 2014a|b ( for com¬ 
plete details. We explore 4 DM models: CDM and WDM 
with mx = 2.25,3 and 5keV. Although we cite values 
assuming thermally decoupled relativistic particles, these 
numbers can also be converted into sterile neutrino masses 
(msteriie^) using ( Viel et al.|2005 l 


"^sterile u 


4.43 keV 


mx y^Yo.l225\ 

IkeV j j’ 


( 1 ) 


yielding msteriieu = (13.1,18.6, 36.8) keV corresponding to 
rux = (2.25, 3, 5) keV respectively. 


2.1 Merger trees and the baryonic 
implementation 

We start by constructing 400 merger trees starting at 
2 = 4 linearly distributed across the halo mass range 
log(MA/M 0 ) = 9 — 13 for the 4 DM models considered. 
The merger-trees use 320 equal redshift steps (A 2 = 0.05) 
between 2 = 20 and 2 = 4 with a mass resolution of 
Mres = 10® Mq using the modihed binary merger tree al¬ 
gorithm with smooth accretion detailed in [Parkinson et ah] 
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1 2008 1 and Benson et al. (20131. We scale the relative abun¬ 


dances of the merger tree roots to match the « = 4 Sheth- 


Tormen halo mass function (HMFs; Sheth & Tormen 19991 
and have verified that these yield HMFs in good agreement 
with the Sheth-Tormen | Sheth fc Tormen|1999 1 HMF at all 

2 . 

We implement the merger trees with baryonic physics 
including star formation, supernova (SN) feedback driven 
ejection of gas, and the merger/accretion/ejection driven 
evolution of the gas and stellar masses. Our model is based 
on the simple premise that the maximum star formation ef- 
hciency of any halo is limited by the SN binding en¬ 

ergy required to unbind/eject the rest of the gas and quench 
further star formation, up to a maximum threshold value 


of /* (see Dayal et al. 2014a|. This model has two 2 and 


mass independent free parameters whose values are selected 
to match the evolving UV LF: the maximum threshold star 
formation efficiency (/*) and the fraction of SN energy that 
goes into unbinding gas (/m). While /„ affects the faint-end 
slope of the UV LF where feedback is most effective, /* de¬ 
termines the normalization at the bright-end where galaxies 
can form stars with the maximum allowed efficiency. 

We implement this simple idea proceeding forward in 
time from the highest merger tree output redshift, 2 = 20 . 
At any 2 step, the initial gas mass, Mg^i{z), in a galaxy 
depends on its merger history: while galaxies that have no 
progenitors are assigned a value Mg^i{z) = {Q,bWm)Mh{z), 
the value is determined both by the gas mass brought in 
by merging progenitors as well as that smoothly-accreted 
from the intergalactic medium (IGM) for galaxies that have 
progenitors. A part of this gas forms new stellar mass, 

M* ( 2 ) with the hnal gas mass left depending on the ratio of 
the (instantaneous) energy provided by exploding SN and 
the potential energy of the halo. We note that at any step, 
the total stellar mass in a galaxy is the sum of mass of the 
newly-formed stars, and that brought in by its progenitors. 

For simplicity, we assume every new stellar population 
has a fixed metallicity of O.OSZq and an age to = 2 Myr. 


Using the popula tion synthesis code STARBURST99 (Lei- 
therCT et al.|l999 l, its initial UV luminosity (at A = 1500 A) 
can be calculated as Lj 7 v( 0 ) = 1 O^^ °^^(M,/M 0 ) ergs“^A“^ 
and the initial output of ionizing photons can be calculated 
as Nion{0) = 10'‘® ®^®®(M*/M©)s“^. Further, the time evo¬ 
lution of these quantities can be expressed as 


Luv{t) = Luv{0) ~ 1.33log-h 0.462 

to 

Nio„{t) = iVion(O) — 3.92 log -I- 0.7 

to 


( 2 ) 

(3) 


For any galaxy along the merger tree, its UV luminosity 
and ionizing photon output rate are the sum of the values 
from the new starburst, and the contribution from older 
populations accounting for the drop with time. We therefore 
have these quantities for 320 2 steps between 2 = 20 and 
2 = 4. 


As shown in Dayal et al. ( 2014a|b I, our model repro¬ 
duces the observed UV LF for all DM models (CDM and 
WDM with nij: = 1.5,3 and 5keV) at 2 ~ 5 — 10 over 2.5 
orders of magnitude in luminosity (7 magnitudes in Muv) 


and predicts the 2 -evolution of the faint end UV LF slope in 
in addition to reproducing observables including the SMD 
and mass-to-light ratios using fiducial parameter values of 
/* = 0.038 and fw = 0.1. We maintain these fiducial param¬ 
eter values in all the calculations carried out in this work. 
Further, we replace 1.5 keV WDM with a heavier mass of 
2.25keV to get tighter constraints on rrix; we note that the 
latter mass naturally fits all of the observed data sets pre¬ 
sented in Dayal et al. ( 2014a|b I. 


2.2 Modelling reionization 


The quantity which plays a crucial role in modelling reion¬ 
ization is the rate of ionizing photons per unit comoving 
volume hion (s~^ Mpc“^) in the IGM, which is essentially 
obtained by (i) integrating the photon production rate Nion. 
of individual galaxies over relevant halo mass ranges, and 
(ii) multiplying by the fraction fesc of ionizing photons that 
escape from halos into the IGM. 

Although reionization is driven by the H i ionizing pho¬ 
tons produced by early galaxies, the ultra-violet background 
(UVB) built up during reionization suppresses the bary¬ 
onic content of galaxies by photo-heating/evaporating gas 
at their outskirts (Klypin et al. 1999| [Moore et al. 1999 


jSomervillej [20021 preventing efficient cooling and slowing 
down its progress. In order to account for the effect of reion¬ 
ization feedback on hion, we assume total photo-evaporation 
of gas from halos below Mmin = 10®Mq (at any 2 ) within 
ionized regions. In other words, halos below Mmin and can 
neither form stars nor contribute any gas in mergers if they 
are forming inside the regions which have already been ion¬ 
ized. Hence the globally averaged hion is given by 


hion{z) = fosc{z)[Qii{z)hii{z) -I- [1 - Qii{z)]hi{z)], (4) 

where Qii{z) is the volume Hlling fraction for ionized re¬ 
gions, and hi I (hi) is the photon production rate density 
within ionized (neutral) regions. Clearly hi contains con¬ 
tribution from all sources which can form stars and is un¬ 
affected by the UVB, while hn represents the case where 
sources below Mmin do not contribute to ionizing photons. 
At the beginning of the reionization process the volume 
Hlled by ionized hydrogen is very small (Qii « 1 ) and 
most galaxies are not affected by the UVB, thus giving 
hion{z) ~ hi(z). However, as Qu increases and reaches a 
value ~ 1 (note reionization is said to be complete when 
Qii = 1), all galaxies less massive than Mmin are feedback- 
suppressed such that hion{z) ~ hii{z). 

The reionization history, expressed through the evolu¬ 
tion of Qii, can then be written as 

dQ 11 diiion Qii dt 

dz dz tree dz ’ 

where dt/dz = [H{z){l + 2 )]“^. Here, the first term on the 
right hand side represents the ionizing photon production 
which drives reionization, and the second term shows the 
decrease in the H ii volume filling fraction due to recombi¬ 
nation of electrons and protons to form H i . Further, nn 
is the comoving hydrogen number density, and tree is the 
recombination time that can be expressed as 
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Figure 1. The CMB electron scattering optical depth (res) as a function of redshift for the 4 DM models considered in this paper, as 
marked in each panel. In each panel, the 4 lines show results using four different values of fesc = 1-0 (dot-dashed line), 0.5 (dashed line), 
0.35 (dotted line) and the fiducial ^-dependent value marked (solid line); we use the same ^-dependence of fesc for both the 3 and 2.25 keV 
cases. The 2 -dependent fesc value has been obtained by simultaneously fitting to Tes and ionizing photon emissivity observations. The 
horizontal dashed line show s the central value for Tes inferred by Planck with the shaded region showing the errors allowed | |Planck| 
Collaboration et al.||2014b l. As seen, while a 2 -dependent fesc gives reasonable results for CDM and mx ^ 3keV WDM, the 2.25 keV 


model is effectively ruled out by current CMB constraints, even using the maximum allowed value of fesc = 1- 


1 

Xnjf (1 + z)^aB C 


( 6 ) 


where as is the hydrogen case-B recombination coefficient, 
X = 1-08 accounts for the excess free electrons arising from 
singly ionized helium and C is the IGM clumping factor 


which we assume to be evolving as (Pawlik et al. 2009 
Haardt fc Madau|2012 ) 


< > 

< riHii 


l + 43z"^'^\ 


(7) 


As seen above, the photon production rate depends on 
fesc and hence we need to know the value of the param¬ 
eter to calculate the reionization history. In this work we 
assume fesc to be independent of the halo mass, and take to 
be a function of only 2 . Once the function fesc{z) is chosen, 
the reionization history of the model can be worked out. As 
shown in Sec. |3.1| that follows, simultaneously fitting to ob¬ 
servations of the CMB electron scattering optical depth (Tea) 
and the ionizing emissivity (hion) require a z-dependent fesc 
that evolves as 


fesc = fo 


l + z 


( 8 ) 


where fo and a are z-independent parameters whose values 
for different DM models are shown in Table [T] 


3 RESULTS 

Now that the model has been described we start by dis¬ 
cussing the model constraints on rrix using observed values 
of the CMB electron scattering optical depth (Tea) and ion¬ 
izing emissivity (hion). We then present the inferred reion¬ 
ization history and sources for all the 4 DM models, as de¬ 
scribed in what follows. 


3.1 Joint reionization constraints from the CMB 
optical depth and UVB emissivity 


Combining Planck and WMAP low-/ polarization data, the 


current best estimate of Tea = O.OSQIq q}^ (Planck Collab- 
oration et al.||2014b|). We calculate Tea at any given z by 


solving the equation: 

Tes{z) = arc / dtne(l + z)^ 
JO 


(9) 


where ne(z) = Qii(z)nH is the global average comoving 
value of the electron number density and ot ~ 6.6524 x 
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Figure 2. Redshift evolution of the H i ionizing photon emissivity (for all theoretical galaxies) for the 4 DM models considered in this 
paper, as marked in each panel. In each panel, the 4 lines show results using four different values of fesc = 1-0 (dot-dashed line), 0.5 
(dashed line), 0.35 (dotted line) and the fiducial 2 :-dependent value (solid line); we use the same ^-dependence of fesc for both the 3 and 
2.25 keV cases. As clearly seen, a constant fesc value severely over-predicts hion compared to the observations (points) at 2 ~ 5,6 in 
all DM models. The observational results (and associated error bars) have been calculated following the approach of |Kuhlen fc Faucher- 


|Giguere l|2012||, i.e., by combining the observational constraints on Fhi from|Wyithe & Bolton|||2011|l with Xmfp from|Songaila & Cowie 
1 |2010[ |. See text for details. 


Model 

fo X 100 

a 

CDM 

1.65 

5.4 

nix = 5 keV (13.1 keV) 

1.25 

6.8 

nix = 3 keV (18.6 keV) 

1.2 

8.2 

rrix = 2.25 keV (36.8 keV) 

1.2 

8.2 


Table 1. The parameter values for 2 -evolution of the escape frac¬ 
tion fesc for different DM models; the numbers in brackets show 
the sterile neutrino mass corresponding to The 2 -dependence 
of fesc is taken to be fesc{z) = /o[(l + ^)/7]“. 


10 is the Thomson scattering cross-section. Note that 

fesc governs the evolution of Qn as shown in Eqns. [^and 

0 


A second observable that needs to be fit by any reion¬ 
ization model is the emissivity of ionizing photons at 2 < 6 . 
It is currently believed that reionization was completed in 
a “photon-starved” manner, which is constrained mainly 
by the observations of H i photoionization rate Fhi from 


quasar absorption lines (Lya forest) at 5 ^ 2 ^ 6 (Bolton 
fc Haehnelt||2007 K In order to calculate Fhi in reionization 


models, one needs to model the mean free path Xmfp of ion¬ 


izing photons in addition to hion (see, e.g., Choudhury & 


Kuhlen & Faucher-Giguere| (|2012 

and combine the obser- 

vational constraints on Fhi from 

Wyithe & Bolton (20111 

with Xmfp from Songaila & Cowie (20101 to obtain the ob- 


from Kuhlen & Faucher-Giguere (20121: 7 = 1 for the source 
spectral index and /3 = 1.3 for the H i column density dis¬ 
tribution. Varying these within allowed would only result in 
larger error bars, hence leaving our results unchanged (see 
Kuhlen fc Faucher-Giguere|2012 1 . 


The most common approach adopted in calculating 
the reionization history (and Tea) is to assume a redshift- 
independent constant value of feac (e.g. Schultz et al.|20l4 |. 
We start by following this basic approach to calculate Tea 
and hion for our 4 DM models, to compare to observations. 
As shown in Fig. matching to the Planck optical depth 
constraints requires feac ~ 0.5 for GDM and mx it 3keV 
WDM models; a value of feac = 0.2 does not result in 
enough electrons to produce the measured optical depth at 
high- 2 . As shown in Dayal et al. 12014b see Fig. 3), struc¬ 
ture formation is progressively delayed with decreasing rux, 
resulting in a lack of H i ionizing photons which naturally 
leads to a low value of Tea- Indeed, our model shows that 
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Figure 3. The redshift evolution of the fiducial value of fesc = 
used in this work (see text for details). The results for 
2.25 and 3keV are the same since we use the same z-dependence 
for both. As seen, progressively lighter WDM particles require a 
larger photon escape fraction at earlier times to compensate for 
the lack of low mass halos. 


Tes is lower than the Planck limits even in the case that 
fesc = 1, where all the H i ionizing photons escape out of 
galaxies and contribute to reionization ruling out WDM with 
rux ^ 2.25 keV0 Although they match the observed Tes, z- 
independent constant values of fesc = 0.5(1) severely over¬ 
predict the ionizing photon emissivity by about 1.6 (1.9) 
orders of magnitude at z ~ 6 as shown in Fig. A sce¬ 
nario wherein /esc is z-independent is therefore ruled out 
by the emissivity constraints, irrespective of the DM model 
considered. 

Reconciling these two data sets requires a z-dependent 
fesc = fQ[{l + z)/7]°‘with a > 0, which provides enough pho¬ 
tons at high-z to obtain the right optical depth Tes, whilst 
yielding reasonably low values at later times; we choose 
/o and a for each DM model such that both the constraints 
are satisfied simultaneously. The fact that one requires ion¬ 
izing sources with higher efficiency to match the two data 
sets has been noted earlier (see e.g., Mitra et al.|201l' 20121. 
With respect to CDM, the number density of the (low mass) 
galaxies occupying the faint end of the UV LF {Muv ~ ~11) 
is increasingly suppressed with redshift as the WDM particle 
mass decreases from 5 to 3 keV. This must be compensated 


^ Note, however, that the constraints on Tes from 

the recent (12/2014) Planck data seem to have been 

(tentatively) revised to Tes = 0.079 ± 0.017, see 

http://www.cosmos.esa.int/web/planck/ferrara2014 for 

details. In that case one may be able to reconcile the nix = 2.25 
keV case with CMB data by choosing fesc ^ 0.3 which yields 
Tes ^ 0.062. 


by an increase in fesc (or a) to obtain the required Tes at 
early epochs: the values of fa and a which best fit the cen¬ 
tral Tes value, while being within hion bounds, for different 
DM models are shown in Table [T] Given that even the max¬ 
imum allowed value of fesc = 1 is not enough to reach the 
measured Tes for the 2.25 keV WDM scenario, we simply use 
the values as in the 3 keV case0 

We show the redshift evolution of the best-fit fesc 
values for all the DM models in Fig. As mentioned, 
with the largest number of low mass halos available to 
provide H i ionizing photons CDM requires a shallower z- 
dependence on fesc- A successive decrease in the number of 
low mass halos from CDM to 5 to 3 keV requires a larger 
fesc at early times, resulting in a steeper z-dependence. To 
quantify, we find fesc — 1.4% at z ~ 6 which rises to 50% 
at z ~ 12,11,10 for CDM, 5 and 3 keV WDM, respectively. 


Reionization histories and sources in different 
DM models 


3.2 


As discussed in Secs. |2.2| and |3.H given the source galaxy 
population, the reionization history (or the global z- 
evolution of Qii) is fixed once the z-evolution of fesc is 
determined. We start by noting that CDM allows galaxy 
formation on the smallest scales, with low mass galaxies be¬ 
ing progressively suppressed as m-x decreases. It naturally 
follows that low-mass galaxies in the CDM model will be the 
most affected by UVB feedback, drawing out the process of 
reionization with this effect decreasing with the WDM parti¬ 
cle mass. This is the behaviour shown in Fig.|^ reionization 
starts earliest in the CDM model and is 50% complete by 
z ~ 13, initially driven by low mass halos (see also Bromm 
fc Yoshida||2011 1. Suppression of star formation in the nu¬ 
merous sources with Mh IO^Mq leads to a “stalling” of 
the reionization process with only a 10% change in Qu over 
the 200 Myrs till z ~ 9. Thereafter reionization extends for 
about 650 Myrs driven by more massive galaxies, with the 
entire volume being reionized by z ~ 5. 

From the same figure we see that while reionization 
starts off at a slower rate for rrix = 5 and 3keV WDM, 
it “catches up” fairly quickly due to the slower z-decline 
of fesc that results in a larger number of ionizing photons; 
indeed, the IGM is 50% by z ~ 11.5 — 12.5 in these mod¬ 
els. As expected, the progressive lack of low-mass galaxies 
means there is a slightly (almost no) stalling of the reion¬ 
ization process in the 5(3) keV WDM scenario resulting in 
a quicker end to reionization: indeed, reionization is over by 
a z ~ 6 (8.5) in the 5(3) keV WDM scenario, having taken 


® It is, however, possible to obtain a Tes 0.062 for the nrix = 
2.25 keVease, while at the same time match the hion constraints 
at z = 6, if we choose fa = 1.2 and a = 8.2. These would then be 
consistent with the latest (12/2014) Planck results. We also find 
that the rox = 1.5 keV WDM gives Tes ~ 0.062 in the unlikely 
scenario that constant fesc = 1, i.e., the new Planck results would 
possibly weaken the lower limit on nix from 2.25 keV to 1.5 keV 
if we assumed every ionizing photon produced by 1.5 keV galaxies 
to have contributed to reionization. 
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Figure 4. Volume filling fraction of ionized hydrogen as a func¬ 
tion of z for the 4 DM models considered using the fiducial model. 
As seen, reionization is progressively delayed (and its rate en¬ 
hanced) for decreasing rrij, when compared to CDM due to pro¬ 
gressive suppression of small scale structure. This naturally re¬ 
sults in reionization being the most “drawn-out” for CDM where 
the effects of the UVB are felt most strongly. 


roughly 780 (420) Myrs, as compared to CDM which re¬ 
quires about 1.02 Gyrs. A delay of about a few tens to a 100 
Myrs in structure formation results in reionization starting 
at a 2 ~ 16.5 in the 2.25 keV WDM scenario. The larger halo 
masses result in this model being independent of the feed¬ 
back effects of reionization, resulting in a smooth 2 -evolution 
of Qii from completely neutral to ionized in the 520 Myrs 
between 2 ~ 16.5 to 6.7. In principle, these different sce¬ 
narios could be distinguished with future 21 cm cosmology 
surveys with next-generation facilities such as the Square 
Kilometre Array (SKA). 

The next question that needs to be answered concerns 
the main reionization sources in these 4 DM models, for 
which we calculate the fractional contribution to the to¬ 
tal ionizing emissivity at 2 ~ 5 from galaxies in varying 
halo mass, stellar mass and UV magnitude bins as shown in 
Fig. HI Starting with the halo mass (panel a), we find that 
CDM galaxies with Mh ~ 1 O®M 0 provide only about 35% of 
the total ionizing photons, with the dominant contribution 
coming from lower halo masses. As noted, small mass halos 
are increasingly suppressed with decreasing in WDM 
models which leads to an increase in the fractional contri¬ 
bution of Mh ~ 1 O®M 0 galaxies. Indeed, their contribution 
increases from 50 to 63 to a 80%, as decreases from 5 to 
3 to 2.25 keV. As expected from the decreasing number den¬ 
sities with increasing Mh, the contribution of Mh Ti, 1 O^®M 0 
halos is lower (~ 4% for CDM) and shows the same trend of 
increasing with decreasing mj,, to ~ 15% for the 2.25 keV 
model; this implies galaxies with Mh — 1 O®“^®'®M 0 pro¬ 


vide about 30% (65%) of the total ionizing photons in CDM 
(2.25keV WDM). From the same panel, we see that the 
UV suppression of low mass halos leads to a flattening in 
their fractional contribution for CDM. However, the frac¬ 
tional contribution of low mass halos rises much more steeply 
with decreasing m^, since halos collapse at scales larger than 
that affected by the UVB. 

We now express the fractional contribution of galaxies 
in terms of observables including the stellar mass, and UV 
magnitude (panels b and c in Fig. [^. We find that galax¬ 
ies with M, 1 O^M 0 contribute 20 % to the total ionizing 
photon emissivity in CDM. Given that forthcoming instru¬ 
ments like the JWST will be able to detect galaxies with 


M* 


1 O''^M 0 (Fig. 4; Dayal et al. 


2014b 


Pawlik et al. 


20111 , our results imply that the bulk (~ 80%) of the stel¬ 


lar mass responsible for producing reionization photons will 
likely not be directly detectable. As explained above, the 
contribution to the total ionizing emissivity rises with de¬ 
creasing rrij;, with M„ Si 1 O^M 0 galaxies contributing about 
~ 65% to the total for the 2.25 keV model. From the same 
panel, we see that galaxies with M* 1 O®M 0 contribute a 
negligible 2 % to the total emissivity, with this value rising 
to ~ 8 % in the 2.25 keV model. 

Finally, we find that galaxies brighter than Muv = —13 
provide ~ 63% of the total ionizing photons for GDM imply¬ 
ing that the rest 37% must come from even fainter galax¬ 
ies. This value rises to about 80% for the 2.25 keV model 
where fewer lower mass/luminosity halos collapse into bound 
structures. Currently detected galaxies {Muv ~ — 18) only 
contribute 6 % of the total emissivity in CDM that rises to 
18% in the 2.25 keV model with the latter showing a steeper 
2 -evolution of the emissivity driven by the faster galaxy as¬ 
sembly. 

To summarize, we find the bulk of the reionization pho¬ 
tons come from galaxies with Mh ~ 1 O®M 0 , M* 1 O^M 0 
and —18 Muv ~ — 13 in CDM. The progressive sup¬ 
pression of low-mass halos with decreasing m^, leads to 
a shift in the “reionization” population to larger (halo 
and stellar) masses of Mh Si 1 O®M 0 and M, IO^Mq for 
nix Sj 3keV WDM, although the UV limits effectively re¬ 
main unchanged. 


4 CONCLUSIONS AND DISCUSSION 

Although the Cold dark matter (CDM) paradigm has been 
enormously successful in explaining the large scale struc¬ 
ture of the Universe, it shows a number of small scale prob¬ 
lems that can be alleviated by invoking (keV) Warm Dark 
Matter (WDM) particles that erase small scale power, sup¬ 
pressing the formation of low-mass structures. Given that 
the low-mass high -2 galaxies are believed to be the main 
contributors to the process of cosmic reionization, our aim 
is to use reionization tracers including the Planck GMB 
electron-scattering optical depth (tes) and the measured ion¬ 
izing photon emissivity (hion) to put constraints on the mass 
{nix) of WDM particles. We use a merger-tree based semi- 
analytic model that traces the DM and baryonic assembly of 
high -2 (2 ~ 5 — 20) galaxies in 4 DM cosmologies: GDM and 


© 0000 RAS, MNRAS 000, 000-000 










Figure 5. Fractional contribution to the cumulative H I ionizing photon density at 2 = 5 by galaxies of different halo masses {Mh\ left 
panel), stellar masses (M»; middle panel) and UV magnitude bins {Muy; right panel) for the fiducial DM models considered in this 
paper, as marked. Left panel: the solid and dashed lines show the contribution from galaxies with x ^ 10® and respectively. 

Middle panel: the solid and dashed lines show the contribution from galaxies with x 10^ and IO^Mq, respectively. Right panel: the 
solid and dashed lines show the contribution from galaxies with x — 13 and —18, respectively. 


WDM with rux = 2.25, 3 and 5keV. Our model includes the 
key baryonic processes of star formation, SN feedback and 
the merger/accretion/ejection driven evolution of gas and 
stellar mass, and has been shown to fit observables over 7 
magnitudes (including the UV LFs, M/L ratios and SMD) 
using only two 2 and mass-independent free parameters with 
fiducial values of 3.8% for the maximum star formation effi¬ 
ciency and 10% of the SN energy going into unbinding gas 
(Dayal et al. |2014a|b| |. In this work, we use the stellar mass 
assembly driven ionizing photon output from these theoret¬ 
ical galaxies to model the reionization history. We include 
the effect of external feedback, i.e., an evolving homogeneous 
UVB suppressing the baryonic content of galaxies below a 
threshold mass {Mh — 10®Mq) to explore the WDM masses 
allowed by reionization constraints. The only additional pa¬ 
rameter introduced is the fraction of ionizing photons (/esc) 
that escape out of galaxies and ionize the IGM. 

We find that ~ 2.25 keV WDM models are ruled out 
because they yield Tes values lower than that measured by 
Planck even in the (unlikely) case that all the H i ionizing 
photons produced by every galaxy escape into the IGM (i.e. 
fesc = !)• Thus, using a model that incorporates the key 
physics driving the baryonic content of high -2 galaxies (star 
formation/feedback/mergers) yields a limit of rux fit 3keV 
that is competitive with r esults of 3.3 keV obtained 

using the Lyman-a forest ( Viel et al.|2013[ )P] 

We also find that while constant values of /esc ~ 0.5 
fit the observed Tes for GDM and 3keV WDM, a 

constant /esc value severely over-predicts the emissivity by 
about 1.6 (1.9)orders of magnitude for /esc = 0.2(0.5). Rec- 


^ The constraints obtained from reionization models could, how¬ 
ever, weaken to ms, 1.5 keV in case the bounds on Tes are re¬ 
vised, e.g., as in the latest (12/2014) Planck results. 


onciling these two data sets requires a z— dependent fesc = 
/o(^^)“ which provides enough photons at early times to 
obtain the right Tes value, whilst yielding low hion values at 
later times. We find that simultaneously fitting to the two 
reionization data sets requires /o = (1.65,1.25,1.2)% and 
a = (5.4, 6 . 8 , 8.2) for GDM and mx = 5,3keV WDM, re¬ 
spectively. We use the same values for 2.25 keV as for 3keV 
since the former never reaches the required Tes value. Quan¬ 
titatively, we find fesc — 1.4% at 2 ~ 6 for all models, rising 
to 50% at 2 ~ 12,11,10 for GDM, 5 and 3keV WDM, re¬ 
spectively. 

As a result of the presence of the most small-scale struc¬ 
ture, reionization starts earliest in CDM, and is 50% com¬ 
plete by 2 ~ 13. Thereafter, the rising UVB suppresses star- 
formation in the very halos that drive reionization drawing 
out its end stages. Indeed, reionization is “photon-starved” 
in GDM and takes a total of about a Gyr to go from a com¬ 
pletely neutral to ionized IGM. While the progressive lack 
of small-scale structure leads to a delay in the start of reion¬ 
ization with decreasing mx, it also transpires that existing 
galaxies are less affected by the UVB leading to an enhanced 
rate in its progress. While the latter effect is dominant for 
nix ~ 3 keV resulting in reionization lasting for about 780 
and 420 Myrs for = 5 and 3keV respectively, reioniza¬ 
tion lasts for about 580 Myrs in the 2.25 keV model, given 
its lack of low-mass halos. 

Finally, we find that the bulk of H i ionizing photons 
come from galaxies with Mh ^ 1 O®M 0 , M, IO^Mq and 
—18 /S Muv ^ — 13 in GDM. The progressive suppression 
of low-mass halos with decreasing nix leads to a shift in the 
“reionization” population to larger (halo and stellar) masses 
of Mh ^ 1 O®M 0 and M. > 1 O'^M 0 for nix ^ 3keV WDM, 
although the UV limits effectively remain unchanged. 

We end with a few caveats. Firstly, we have ignored 
SN radiative losses that could significantly reduce the to- 
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tal energy available to drive winds. However, a decrease in 
this total energy could be countered by scaling up the frac¬ 
tion of the total energy we put into driving winds from the 
fiducial value of 10% used in this work. Secondly, |Pawlik| 


et al. (20091 have shown that UVB photo-heating reduces 


the clumping factor of the IGM since the additional pressure 
support from reionization smoothes out small-scales density 
fluctuations. While we have used their results for an over¬ 
density of 100, we have confirmed that our results do not 
change using threshold values of 50, or 200. Thirdly, while 
we assume that all halos with mass Mh ^ 10®Mq are feed¬ 
back suppressed, we find our results are equally consistent 
with observations (ves and hion) using values ranging be¬ 
tween 10*'® — 1O®'®M0. Finally, we have made the simpli¬ 
fying assumption of using a mass-independent f^sc at all z. 
This is partly motivated by the uncertainty regarding the 
mass-dependence of /esc: while some authors find /esc to 


decrease with an increase in the halo mass (Razoumov & 
Sommer-Larsen||2010| |Yajima et al.|[^011| [Ferrara fc Loeb' 


2013), other works have found the opposite trend (Gnedin 
et al.|2008 Wise fc Cen|2009 1. Although there are a number 


of caveats involved, our model hts observed UV LFs at faint 
end validating the theoretical underlying stellar population. 
Unless the results on LEG faint-end or Tcs evolve signifi¬ 
cantly, our results rule out rrix ^ 2.25 keV WDM models. 
Significant progress is expected to be made by comparing 
our model predictions with actual faint LEG data such as 
that imminently expected from the Frontier Fields, and from 
forthcoming observatories such as the JWST. Furthermore, 
if the claim of 3.5 KeV X-ray line holds to be true, it will 
usher in an exciting era for WDM dominated cosmology. 
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